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The Ge-C local mode in epitaxial GeC and Ge-rich GeSiC alloys 

W. H. Weber** 

Physics Department, MD-3028/SRL Ford Motor Company, Dearborn, Michigan 48121-2053 

B.-K. Yang and M. Krishnamurthy 

Department of Metallurgical and Materials Engineering, Michigan Technological University , 
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(Received 19 January 1998; accepted for publication 27 May 1998) 

The Raman signature of the local Ge-C mode for substitutional C is identified as a narrow line (8 
cm” 1 full width at half maximum) near 530 cm” 1 in alloy films of Gej _ V C V (0^y^0.07) grown 
on Ge (100) substrates by low-temperature (200 °C) molecular beam epitaxy. The intensity of :he 
Ge-C line relative to the c-Ge line suggests that only a small fraction of the nominal C is in 
substitutional sites. In ternary alloys of Ge| _ x _ v Si x C v with jc = 0.1 and 0.2 and y = 0.03, the Ge-'C 
mode disappears, suggesting a strong bias towards C bonding with Si as opposed to Ge. In 


Ge|_ x Sn x films the Ge-Sn mode is seen at 263 cm 
[S0003-6951(98)00231-9] 

The addition of small amounts of C to the Ge-Si alloy 
system is considered an attractive method for reducing the 
lattice-constant mismatch in layered growth of band gap en¬ 
gineered materials. 1 The introduction of C also offers an ex¬ 
tra degree of freedom to tune the electronic properties by 
varying both strain and composition. Although the solubility 
of substitutional C in bulk Si or Ge is extremely low, 2 the 
solubility near a surface can be several orders of magnitude 
higher, 3 which offers the possibility of producing metastable 
films with super-saturated concentrations as high as 2-3 
at. 9c 4 Thin-film samples with concentrations this high have 
been achieved by a variety of techniques including molecular 
beam epitaxy (MBE) 5-8 and chemical vapor deposition. 9 " 11 

In Si and Si-rich alloys, substitutional C produces a char¬ 
acteristic Si-C local Raman mode at 605 cm" 1 . 5,6,12 ' 15 Ob¬ 
servations of this mode are significant, since they confirm the 
presence of substitutional C, as opposed to interstitial or sur¬ 
face C, and they yield microscopic information regarding the 
distribution of C in the lattice. 12 ’ 15 " 18 The mode intensity can 
be used to estimate the absolute concentration of substitu¬ 
tional C, and the satellite structure around the main peak can 
indicate the occurrence of preferential instead of random sec¬ 
ond or third nearest-neighbor C-C spacings. 

In this letter we report on the corresponding Ge-C local 
mode in Ge and Ge-rich alloys. Although the existence of 
this mode has been anticipated in several recent 
papers, 61319-20 and its position even tentatively assigned, it 
has not been unequivocally identified due to its extreme 
weakness in Si-rich alloys or because of interference with 
other Raman features from the film or the substrate. As with 
the Si-C local mode, observations of the Ge-C mode con¬ 
firm the presence of substitutional C and shed light on the 
microscopic C distribution. The latter feature is particularly 
important in the Ge-Si-C alloys, since Kelires has shown 
that Ge-C bonds are energetically unfavorable compared 
with Si-C or Ge-Si bonds, so that highly nonrandom distri¬ 
butions can be expected. 21 

The samples in this study are epitaxial films of 
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Gej. x . y Si x C y alloys, nominally 50-70 nm thick, grown on 
Ge (100) substrates using a Riber 32 MBE system with base 
pressure lower than 5X10” 10 mbar as described in more 
detail elsewhere. 20 The films were analyzed in situ by reflec¬ 
tion high-energy electron diffraction (RHEED) and ex situ 
by x-ray diffraction (XRD), transmission electron micros¬ 
copy (TEM), and atomic force microscopy (AFM). Ruther¬ 
ford backscattering (RBS) was used to calibrate the various 
Ge, Si, and C sources. The Ge(100) substrates were de¬ 
greased and cleaned in de-ionized water and then In-bonded 
to a Mo block. The oxide was desorbed in situ at 500- 
550 °C and a 20-nm Ge buffer layer was grown at — 1 nm/ 
min and 200 °C substrate temperature. A sharp (2X1) 
RHEED pattern from clean Ge was observed in all samples. 
The Si was deposited either from a Knudsen source with a 
TaC crucible or from a rod-fed electron beam source. The C 
was deposited from a rod-fed electron beam source, whose 
flux was varied to control the C conc^Tjratior. The ratio of 
this flux to that of the Ge flux we refer to as the “nominal” 
C concentration. 

During codeposition of the Ge and C (and Si, in the case 
of SiGeC growth), the RHEED pattern undergoes significant 
changes. For all samples the (2X1) pattern transforms gradu¬ 
ally to a (1 Xl) pattern followed by three-dimensional island 
formation and strong {311} faceting during later stages of 
growth. The higher the C concentration, the sooner these 
changes occur. Plan-view and cross-sectional TEM images 
of the 5%C Ge-C film have already been reported. 20 The 
plan-view images show strain contrast with no indications of 
dislocations. Cross-sectional views show {311} facets deco¬ 
rated with C clusters. AFM images of the surface show uni¬ 
formly distributed islands of average size — 50 nm and 
heights of 10 nm or less, which is consistent with the TEM 
results. Attempts to determine the film’s lattice parameter by 
26 XRD scans showed no clear evidence of a peak from the 
epilayer separate from The strong Ge substrate signal. 

Raman spectra were obtained using a Renishaw 1000 
Raman microscope wiui 633 nm excitation. At this wave¬ 
length the films are sufficiently transparent that their Raman 
signals are comparable to those from the substrate. Since the 
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FIG. 1. Raman spccirum of a 50 nm thick, GeC alloy film with nominal 2% 
C concentration. Experimental conditions: 3 mW and 20 min exposure. 

Ge-C local mode occurs near the crystalline Si (c-Si) line at 
521 cm’ 1 , it is crucial that we use a Ge substrate. In all of 
our spectra taken earlier on Si substrates, the c-Si line com¬ 
pletely obscures the much weaker Ge-C local mode. 20 

Figure 1 shows the Raman spectrum from a 50 nm thick 
Ge, _ s C y epilayer with a nominal C concentration of 2%. The 
spectrum is recorded with z(xy)z scattering geometry, for 
which the first-order c-Ge line and the Ge-C local mode are 
allowed and the second-order scattering is largely sup¬ 
pressed. In the magnified plot, the broad featu e peaking near 
575 cm" 1 is the second-order band from die Ge optical 
branch. All the features in this spectrum, excepting the nar¬ 
row line at 530 cm’ 1 , indicated in the figure, also appear in 
the spectrum from a pure Ge (100) wafer. We assign this 
extra line to the Ge-C local mode. This line is absent from 
the Ge-wafer spectrum, and its strength increases with C 
content. It follows the same selection rules as the c-Ge peak, 
i.e.. it is observed in xy and x'x' scattering and is quenched 
in xx and r'v' scattering. In addition, its frequency is con¬ 
sistent with that expected from the known position of the 
Si-C local mode at 605 cm -1 . Compared with the Si-C 
mode, the larger reduced mass of the Ge-C pair will shift the 
frequency down to —546 cm" 1 , which corresponds to a fac¬ 
tor (m Sl _ c /m Ge _ c ) ]/2 , where the w,-./s are reduced masses 
(m~J =m~ ] + m~'). The further downward shift of —3% 
is likely caused by the larger lattice constant for Ge com¬ 
pared with Si, which increases the Ge-C bond length and 
decreases the force constant. 

The Ge-C mode in Fig. 1 has a full width at half maxi¬ 
mum of 8 cm" 1 . The narrowness of this mode is consistent 
with its frequency—well above the top of the one-phonon 
continuum and just below the overtone band from the c-Ge 
optical branch, in a region where the second-order scattering 
is weak. 22 Thus, one-phonon decay processes are forbidden 
and two-phonon processes should be very weak. 



FIG. 2. Spectra of films with different C concentrations and a pure Ge wafer 
recorded under identical conditions as those for Fig. 1. Successive spectra 
are shifted up by 200 counts for clarity. The inset shows the integrated 
intensity of the Ge-C local mode vs nominal %C. The dashed line is a fit to 
a( 1 — exp{— bx}) drawn as a guide to the eye. 

Figure 2 shows an expanded plot of the Ge-C line for 
samples with different C content along with a Ge-wafer 
spectrum. The linewidth increases about 25% at the highest 
C concentration, but the integrated line strength, shown in 
the inset, increases by only about 50%. This result suggests 
that the substitutional C content saturates at a value well 
below the concentration expected from the growth process. 

There are also subtle changes in the c-Ge parts of me 
spectrum as well. The first-order peak broadens slightly 
( — 10%) and decreases in integrated intensity ( — 20%), and 
the structure on the edges and near the peaks of the second- 
order features tends to become less sharp. These effects are 
likely caused by the faceting and reduction in long-range 
order inherent in the growth process. Unfortunately, these 
changes make it impossible to form difference spectra be¬ 
tween C-doped and pure Ge samples without introducing 
extra structure throughout the region below 600 cm" 1 . 

We observe a broad band, peaking near 1400 cm -1 , 
that increases with C content. This band is the signature for 
amorphous C (a-C), 20 ’ 23 indicating that at least some of the 
extra C goes into o-C. There are also additional weak, sharp 
lines in the spectra (near 950 and 1050 cm -1 ) of unknown 
origin, which may be related to specific C or Ge-C clusters. 

In GeSiC alloy films with 10% and 20% Si and nomi¬ 
nally 3% C, we find the usual Ge-Ge and Ge-Si lines, but 
we find no evidence for the Ge-C local mode. A likely ex¬ 
planation for this effect is that the strong preference for C to 
form Si-C bonds leads to every C being paired with a Si 
atom. 21 

In Sii__ v C v alloys Melendez-Lira el al, ]9 found that the 
ratio of integrated intensities of the Si-C mode at 605 cm - 1 
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RC. 3. Raman spectrum of a GeSn epilayer with nominal 5% Sn, recorded 
under the same conditions as for Fig. 1. 

to the r-Si line at 521 cm" 1 was given by —3.7>\ where v is 
the substitutional C fraction. A theoretical estimate of this 
ratio yielded I.9v, which is within a factor of 2 of the ex¬ 
perimental result. If we apply the same theoretical arguments 
they used to estimate the corresponding ratio between the 
Ge-C and r-Ge modes, we predict a ratio -~2.6v for our 
samples. Our measured ratio for the nominal 2 % C sample of 
Fig. 1 is 8 X I0' 4 . which is nearly two orders of magnitude 
smaller than the ^dieted value of 5 X 10 -2 . We conclude 
either that there a strong resonant enhancement of the c-Ge 
mode relative to the Ge-C mode or the substitutional C con¬ 
centration is extremely small. A resonance enhancement can¬ 
not be ruled out. since the 633 nm laser excitation is near the 
E, critical point of the r-Ge electronic response. However, it 
is worth noting that the C concentration suggested by the 
Raman intensity ratio. ~3X 10" 4 , would imply (using Veg- 
ard's law) a tensile strain of - 1 X 10“ 4 . which would be too 
small by about a factor of 10 to be detectable in our x-ray 
studies. 

As a check on the validity of using Raman intensity 
ratios to estimate dopant concentrations in Ge. we examined 
a 50 nm thick Ge-Sn alloy film grown in the same fashion as 
the Ge-C films, in a Ge lattice it is generally easier to get Sn 
into a substitutional site than C. 24,25 The Raman spectrum for 
a nominal 5<7r Sn alloy is shown in Fig. 3. The Ge-Sn vi¬ 
bration is identified as the peak at 263 cm” 1 , and, as ex¬ 
pected, there is no suggestion of a peak in the 450-550 
cm" 1 region. The integrated intensity ratio of this line to the 
r-Ge line is 5 X 10" 2 , which is the same order of magnitude 
as die' nominal Sn fraction, and which is consistent with the 
results on the Si-C mode in Si,_ v C v alloys. 19 This result 
supports the conclusion that our substitutional C concentra¬ 
tion is well below the nominal C fraction. 

The Ge-Sn mode is not truly ’’local” in the same sense 


as the Ge-C mode. Since Sn is heavier than Ge. the Ge-Sn 
vibrations will strongly mix with the continuum from the Ge 
optical branch, resulting in a much broader line. As evidence 
of this mixing, the relative linewidth (Aa>/o>) for the Ge-Sn 
mode is ~3x larger than that for the Ge-C mode. 

In summary , we have identified the Ge-C local mode in 
Ge, -yC y epitaxial films on Ge (100) substrates as a sharp 
line at 530 cm" 1 . The intensity of this line relative to the 
c-Ge line and its variation with C content suggest that the 
substitutional C concentration is well below the nominal C 
fraction. At least some of the extra C appears in the Raman 
spectra as a- C. In GeSiC alloys the Ge-C line disappears, 
indicating a strong preference for the formation of Si-C 
bonds as opposed to Ge-C bonds. 

The authors wish to thank Barry Wilkens (Arizona State 
University) for RBS measurements. Financial support for 
part of this research (M.K. and B.Y.) has been provided by 
the Office of Naval Research (NOOO14-96-1-0793). 
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Microstructural development and optical properties of epitaxial 
Gei _ X C X alloys on Si(100) 

Houghton, Michigan 49931 
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constants show small systematic chang features attributable to Ge-C vibrational 
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modes could be identified using Raman spectroscopy. ^ 
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Silicon-germanium alloys have been studied in order to 
exploit the promise of new devices obtainable by “band-gap 
engineering” on silicon. 1 Recently, there de¬ 

creasing interest in SiGe alloys containing carbon. One ot 
the reasons for substitutionally alloying with C is that (as¬ 
suming Vegard’s law) approximately a 9:1 Ge:C ratio wil 
lattice match SiGeC alloys to Si, circumventing critical 
thickness limitations inherent in SiGe technology. However, 
the relative insolubility of C in both Si and Ge and the pre¬ 
cipitation of various forms of SiC at higher grow 
temperatures 1 are obstacles to the fabrication of these alloys 
Despite these problems, low concentrations of substitutional 
C in the alloys have been successfully incorporated using a 
variety of techniques including chemical-va^r deposition 
(CVD) 4 5 and molecular-beam epitaxy (MBE). 

The Ge-C system is an alternative pathway for Si-based 
band-gap engineering. In addition to the fact that it is a sim¬ 
pler binary system, problems associated with compound for¬ 
mation can be avoided since no stable GeC phase exists^ 
Several researchers have reported oif the growth of 
Ge C, alloys’* " on Si. Previous MBE studies of 
Ge I C,7 1 films were typically performed at higher tempera¬ 
tures (600 °C) on Ge buffers 8 or using Sb surfactants, with 
little emphasis on microstructural development or the dielec¬ 
tric function. Our earlier study of Ge, . X C X e P ltax y on Sl had 
shown that high-C films (>20%) grown at room temperature 
and annealed at higher temperatures formed 3D islands and 
the Ge and C were phase separated. 

In this letter, we report on MBE growth of Ge,-,!-* 
alloys on Si(100) at 200 °C, for C concentrations between 
0% and 10%. The films were characterized in situ by reflec¬ 
tion high-energy electron diffraction (RHEED) and ex situ 
by transmission electron microscopy (TEM), x-ray diffrac¬ 
tion. Raman spectroscopy, and spectroscopic ellipsometry 

(SE). 

The Ge, . X C X films were grown in a Riber 32 MBE sys¬ 
tem with base pressure better than 5 X 10' mbar. Si(100) 
substrates were chemically cleaned and oxidized by boiling 

•’Electronic mail: mohan@mtu.edu 
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in a H 2 0 2 :H->S0 4 mixture. After the sample had been de¬ 
gassed in the MBE chamber at -200 °C overnight, theoxide 
was desorbed by heating to just below 900 °C for about 15 
min Sharp (2X1) diffraction spots indicative of a clean sur¬ 
face were typically observed in RHEED (in some cases a 
couple of monolayers of Ge were deposited on the oxide 
before desorption to produce a better quality surface). Ger¬ 
manium was deposited from a previously calibrated pyrolyttc 
boron nitride Knudsen source at -10 A/min. The C was 
deposited from a rod-fed electron-beam source, and the flux 
waV varied to obtain various C concentrations. Rutherford 
backscattering spectroscopy (RBS) (including C resonance 
was used to measure the C concentrations. An experiment 
consisted of the following steps. After the sample had stabi¬ 
lized at -200 °C (±25 °C), the Ge shutter was first opened 
followed by the C shutter after a 10-30 s delay. The nomma 

thickness of most films was—50 nm. • 

During deposition of the Ge-C alloy, the RHEED p 
tern undergoes changes depending on the C content^ as 
shown in Figs. 1(a)-1(d). At low C concentrations (0%- 
2%), the (2X1) reconstruction is stable and the streakiness is 
indicative of a 2D layer growth with an atomically rough 
surface [Figs. 1(a) and 1(b)]. For medium C films (3%-5%), 
the initial 8-10 monolayers (ML); (1 ML-1.4 A) shows a 
streaky (2X1) pattern, followed by evidence of transmission 
spots indicative of 3D island growth [Fig. 1(c)]. For higher C 
films extra spots indicative of twinning are also apparent in 
the RHEED pattern [Fig. 1(d)], and the films tend to become 
polycrystalline at the later stages of growth. 

Figures 2(a) and 2(b) are cross-sectional TEM images of 
Ge-C films with low C concentration (l%-2% C). Figure 
2 ( a ) (~i% C film) shows that the film is of good epitaxial 
quality with the majority of the misfit defects localized near 
the interface. The high-resolution image in Fig. 2(b) ( i*> 
C) shows the Ge-C/Si interface with misfit defects, a few 
planar defects (twins, stacking faults), and localized regions 
of the lattice where the fringes appear disordered. Figures 
3(a) and 3(b) show plan-view and high-resolution cross- 
sectional TEM images, respectively, of a 50-nm-thick, Ge-C 
alloy with -5% C. The plan-view image [Fig. 3(a)] indicates 
that the film is relaxed with short interfacial misfit disloca- 
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FIG. I. Typical RHEED patterns observed for (a) pure Ge film, (b) film 
with <2Cr C. (c) film with ~5 % C showing transmission spots, and (d) film 
with -109c C showing twinning spots. 


tions as evidenced by the presence of moire fringes and from 
weak-beam dark-field images (not shown). The cross- 
sectional image shows faceting that is consistent with the 
RHEED observation of 3D islanding. There is no clear evi¬ 
dence for the formation of C clusters or precipitates. 

X-ray 20 scans were obtained in the (400) range of Si 
and Ge. The ratio of rW^aoosi is used as a measure of the 
lattice parameter of the relaxed film, and is plotted as a func¬ 
tion of C concentration in Fig. 4. A marginal decrease in the 
Ge lattice parameter with increasing nominal C concentra¬ 
tion is observed. It is estimated (based on Vegard’s law) that 
on an average, a maximum of ~\% C may have been sub- 



FIG. 2. High-resolution cross-sectional TEM images of Ge,-,*:* films 
taken along a <110) zone for (a) film with nominally -19c C and (b) film 
jwith -29c C. Arrows indicate the position of the Ge-C/Si interface. 


FIG. 3. (a) Plan-view and (b) cross-sectional TEM images of a 5% C, Ge-C 
film on Si. The plan-view image was taken under two-beam bright-field 
conditions and shows uniform moire fringes. 

stitutionally incorporated in the films under the growth con¬ 
ditions studied. 

Spectroscopic ellipsometry (SE) measurements were 
performed on the Ge,.^ samples using a variable-angle 
instrument from J. A. Woollam Co. Data were obtained at 
five equally spaced angles of incidence from 60° to 80° and 
over a spectral range of 250—1000 nm. The data were ana¬ 
lyzed to extract the optical constants of the films using a 
three-layer model: Si substrate, Ge^C, layer, and Ge0 2 
layer. The free parameters were the optical constants and 
thicknesses of the Ge | - film and the thicknesses of the 

GeO-. layer. The Ge0 3 index was determined to be ~1.8 
from a separate analysis on pure Ge wafers, and the optical 
comparts of Si were taken from Palik. The free parameters 
wc.e first fit in the long-wavelength region (\>680 nm) 
where there is no sharp structure, by modeling the film opti¬ 
cal constants with a Cauchy formula for the real part of the 
index n and an Urbach tail for the absorption coefficient k. 
The thicknesses resulting from the fit were then constrained 
and the optical constants of the film were found for the full 
wavelength region. This procedure was developed using pure 



FIG. 4. The ratio of the (400) interplanar spacing of Ge to Si (doeJd s ) 
plotted as a function of nominal C content. The arrow indicates the ratio for 
an ideally relaxed pure Ge film on Si. Ratios lower than this value corre¬ 
spond to a reduction in the average lattice parameter of the film due to C. 
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Ge epitaxial layers (with 40-50 A of oxide) for which it was 
found to give good agreement with the results of Aspnes and 
Studna 14 for bulk Ge. 

The most pronounced systematic effect of the C alloying 
is a reduction in the strength of the feature near the E j criti¬ 
cal point (—2 eV), which is shown in Fig. 5. The solid lines 
give the real and imaginary dielectric constants for pure Ge. 
The circles, boxes, and triangles are the data from films with 
nominal C concentration —1%, 2%, and 5%, respectively. In 
addition, less systematic shifts and broadenings also occur at 
higher energies near the E 2 critical point. 

Raman scattering studies were performed on the Ge-C 
films using 633 and 785 nm lasers. All spectra show the first- 
and second-order structure from crystalline Si and Ge. No 
new features attributable to Ge—C bonds could be identified. 
In the samples with 5% C, the Ge first-order line near 300 
cm' 1 showed some broadening, which may result from the 
disorder. In fact, the changes in the electronic structure 
shown in Fig. 5 may result indirectly from an increase in the 
disorder or local strain in the Ge lattice caused by the C, 
rather than from a direct interaction with the C orbitals. 

There are several interesting aspects of microstructural 
evolution of the Ge,.^ films. It is well understood that 
pure Ge grown on Si forms islands (after initial layer 
growth) in order to reduce its strain energy at the cost of 
increased surface area. By maintaining a sufficiently low 
substrate temperature (e.g., 200 °C), the reduced surface dif¬ 
fusion of Ge suppresses island formation? 5 as confirmed by 
Fig. 1(a). 

However, it is seen in our study that the presence of 
C>2%-3% leads to 3D island formation [Fig. 1(c)] even at 
this low substrate temperature. There are several possible 
reasons for the increased islanding tendency in Ge-C films. 
The most obvious one is an increase in the substrate tem¬ 
perature due to radiative heating from the C source. This is 
unlikely because the thermocouple reading from the sub¬ 
strate registers no significant change during the course of the 
experiment. Other possibilities include changes in the sur¬ 
face energy of Ge:C, enhancement of surface diffusion kinet¬ 
ics or a combination of both. While further studies are 
needed to clarify the mechanism, we observe that the growth 
of Ge-C on a relaxed Ge film at higher substrate tempera¬ 
tures also leads to islanding. We therefore suggest that the 
increased driving force for island formation in Ge-C films at 
200 °C may be due to changes in surface energy (i.e., a non¬ 
wetting behavior) induced by the presence of C. 

A second microstructural feature associated with in¬ 
creased C concentration is the appearance of twinning in the 
islands [Fig. 1(d)]. The mechanism for twinning appears to 
be related to the evolution of islands on the surface and may 
be due to the misfits associated with island coalescence. 16 It 
is interesting to note that the microstructural development in 
Ge-C films (i.e., layer, island, twinned island) is qualita¬ 
tively similar to the epitaxial growth of another “immis¬ 
cible" system: Ge-Sn on InSb. 17 

In summary, we show that on Si at 200 °C. epitaxial 
Ge^ x C x alloys with C<29r-3% grow in 2D layers, while 
films with 05 7c form 3D islands after initial layer growth. 
X-ray diffraction indicates that less than \% C may be sub- 
stitutionally incorporated in the films. Spectroscopic ellip- 



FIG. 5. Plot of the real Re(e) and the imaginary lm(e) parts of the dielectric 
function of Ge,..^ as a function of energy for 1% C (circles), 29c C 
(squares), and 5% C (triangles) compared with pure Ge (solid line). Observe 
a systematic change with C concentration near eV (Ej critical point in 
Ge). 

sometry shows a C-related change in the optical constants 
near 2 eV (the E x critical point in Ge) while Raman spec¬ 
troscopy shows no Ge-C-related vibrational modes. 
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(N00014-96-1 -0793). 


‘See for example. J. C. Bean, Proc. IEEE 80. 571 (1992); T. P. Pearsall, 
Crit. Rev. Solid State Mater. Sci. 15. 551 (1989). 

2 For a review see R. A. Soref. Proc. IFEE 81. 1687 (1993). 

3 A. R. Powell, F. K. LeGoues, an^ 3. S. Iyer, Appl. Phys. Lett. 64, 324 
(1994). 

4 Z. Atzmon. A. E. Bair. E. J. Jaqatz, J. W. M.yer, D. Chandrasekhar, D. J. 
Smith. R. L. Hervig, and McD. Robinson. Appl. Phys. Lett. 65, 2559 
(1994). 

5 J. L. Regolini, S. Bodnar. J. C. Oberlin. F. Ferrieu, M. Gauneau, B. Lam¬ 
bert, and P. Boucaud, J. Vac. Sci. Technol. A 12, 1015 (1994). 

6 S. S. Iyer, K. Eberl, M. S. Goorsky, F. K. LeGoues, J. C. Tsang, and F. 
Cardone, Appl. Phys. Lett. 60. 356 (1992). 

7 K. Eberl. S. S. Iyer, S. Zollner, J. C. Tsang. and F. K. LeGoues, Appl. 
Phys. Lett. 60. 3033 (1992). 

8 J. Kolodzey, P. A. Oneill, S. Zhang, B. A. Omer, K. Roe, K. M. Unruh, C. 
P. Swann. M. M. Waite, and S. I. Shah, Appl. Phys. Lett. 67, 1865 (1995). 
9 H. J. Osten, E. Bugiel, and P. Zaumseil, J. Cryst. Growth 142, 322 (1994). 
10 M. Todd, J. Kouvetakis, and D. J. Smith, Appl. Phys. Lett. 68, 2407 
(1996). 

11 M. Krishnamurthy, J. S. Drucker, and A. Challa, J. Appl. Phys. 78, 7070 
(1995). 

12 Handbook of Optical Constants, edited by E. D. Palik (Academic, New 
York, 1985). 

13 See. for example. Guide to Using WVASE32, (J. A. Woolam Co.. 1995). 
14 D. E. Aspnes and A. A. Studna. Phys. Rev. B 27. 985 (1983). 

‘'D. J. Eaglesham and M. Cerullo, Appl. Phys. Lett. 58. 2276 (1991). 
lh See. for example. M. J. Stowell. in Epitaxial Growth. Part B. edited by J. 

W. Matthews (Academic, New York, 1975). 

17 E. A. Fitzgerald, P. E. Freeland, M. T. Asom, W. P. Lowe, R. A. Machar- 
rie, B. E. Weir, A. R. Kortan, F. A. Thiel, Y.-H. Xie, A. M. Sargent, S. L. 
Cooper. G. A. Thomas, anu L. C. Kimmerling, J. Electron. Mater. 20. 489 
(1991). 


2574 Appl. Phys. Lett., Vol. 69. No. 17, 21 October 1996 


Krishnamurthy, Yang, and Weber 




Incorporation and stability of carbon during low-temperature epitaxial 
growth of Ge w C x (x<0.1) alloys on Si(100): Microstructural 
and Raman studies 
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Low-temperature ( — 200 °C) molecular beam epitaxy of Gej_ x C x alloys grown on Si(100) have 
been extensively investigated by in situ reflection high-energy electron diffraction, ex situ x-ray 
diffraction, transmission electron microscopy, and Raman spectroscopy. Carbon concentrations 
were nominally varied from 0 up to —10 at. %. Selected samples were annealed in an Ar ambient 
at 750 °C to evaluate the stability of the thin films. A few films were also grown on Ge substrates. 

With increasing C concentration, the epitaxial growth mode changes from two dimensional layer 
growth to three dimensional island growth. The surface has a tendency to facet along {311} planes 
under certain growth conditions. The microstructure shows an increase in planar defect density with 
increasing C concentration. The x-ray diffraction data show that the lattice parameter decreases with 
increasing C concentration and that a maximum of 1 at. % C is incorporated substitutionally in Ge. 

Raman spectroscopy shows no clear Ge-C signal though extra intensity is measured at the energies 
where Ge-C modes may be expected. Films with nominal C concentrations greater than 2 to 3 at. % 
show clear evidence for amorphous C. We propose that under our growth conditions, nominal C in 
excess of about 2 to 3 at. % remains on the surface as amorphous C and plays an important role in 
3D islanding, defect formation, and {311} faceting during epitaxial growth. © 1997 American 
Institute of Physics. [S0021-8979(97)03719-5] 


I. INTRODUCTION 

Silicon-germanium alloys have been widely studied in 
order to explore the promise of band gap engineering on 
silicon. 1 More recently, there has been a strong interest in 
SiGe alloys containing C. 2 The main reason is that the addi¬ 
tion of suitable amounts of C (a ratio 9:1 Ge:C assuming 
Vegard’s law) is expected to lattice match GeSiC to Si, cir¬ 
cumventing the critical thickness limitations inherent in SiGe 
technology. Furthermore, it is possible to use both alloy con¬ 
centration and strain as a variable for band gap engineering. 
However, there are problems associated with the fabrication 
of these alloys. These include: (i) the relative insolubility of 
C in both Si and Ge and (ii) the precipitation of various 
forms of SiC at higher growth temperatures. 3 Despite these 
problems, low concentrations of substitutional C in the alloys 
have been successfully incorporated using a variety of 
techniques. 4-7 Current research efforts focus on techniques 
to increase the amount of substitutional C in the alloys. For 
example, it has been reported that changes in growth condi¬ 
tions (e.g., growth rates) may change the substitutional-to- 
interstitial carbon ratio in Si-C alloys. 8,9 Most of the epitax¬ 
ial growth studies have been limited to Si-C or Si-Ge-C 
systems. 10 

The Ge-C system is an alternative pathway for Si- based 
band gap engineering. In addition to the fact that it forms a 
model binary system, problems associated with compound 
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formation can be avoided since no stable GeC phase is 
known to exist. Low concentrations of C have been success¬ 
fully incorporated using molecular beam epitaxy (MBE) 11,12 
and chemical vapor deposition. 13 Previous MBE studies of 
Gti- x C x alloys were typically performed at higher tempera¬ 
tures (400-600 °C) on Ge buffers 12 or using Sb surfactants, 11 
with little emphasis on understanding the microstructural de¬ 
velopment or the growth mechanisms in this “immiscible” 
system. Furthermore, the nominal C content of the alloys in 
these studies was limited to x<0.03. 

We have recently reported on the microstructural evolu¬ 
tion and optical properties of Ge } _ x C x alloys (0<x<0.1) 
grown on Si at 200 °C. 14,15 Here we report on a systematic 
study of Ge-C alloys grown on Si as well as on Ge sub¬ 
strates, and the effect of annealing at higher temperatures. 
Microstructural characterization using electron and x-ray dif¬ 
fraction and electron microscopy is correlated with optical 
characterization using Raman spectroscopy. The incorpora¬ 
tion and stability of C in the films and the influence of C on 
the epitaxial growth mechanisms are discussed. 

II. EXPERIMENTS 

The Gej.jC* (0<*<10at. %) films were grown in a 
Riber 32 MBE system with base pressure better than 5 
X 10“ 10 Torn The 2 in. Si(100) wafers were chemically 
cleaned and oxidized by boiling in a H 2 02 :H 2 S 04 mixture 
followed by rinsing in deionized water before being loaded 
into the MBE chamber. After the sample had been degassed 
in the MBE chamber at —200 °C overnight, the oxide was 
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desorbed by heating to just below 900 °C for about 15 min. 
Sharp (2X1) diffraction spots indicative of a clean surface 
were typically observed in reflection high-energy electron 
diffraction (RHEED). In some cases, a couple of monolayers 
of Ge were deposited on the oxide before desorption to pro¬ 
duce a better quality surface. Germanium was deposited 
from a previously calibrated pyrolytic boron nitride Knudsen 
source at — 1 nm/min. The C was deposited from a rod-fed 
electron-beam source, and the flux was varied to obtain vari¬ 
ous C concentrations up to 10 at. %. Rutherford backscatter- 
ing spectroscopy (including C resonance) was used to esti¬ 
mate the C concentrations (with error less than —0.5 at. %). 
A typical experiment consisted of the following steps. After 
the sample had stabilized at —200 °C, the Ge shutter was 
first opened followed by the C shutter after a 10-30 s delay. 
The nominal thickness of Gej_ x C x films was —50 nm. 

In addition to the regular samples discussed above, three 
samples (called A, B, and C) were fabricated under similar 
growth conditions (i.e., 200 °C substrate temperature). 
Sample A had a — 50-nm-thick Ge-C (5% C) thin film 
grown on Si(100) and subsequently capped with —50 nm 
Ge. Sample B was a Ge-C (5% C) film grown on Si(100) 
and capped with — 20-nm-thick Si. Sample C was a —50- 
nm-thick Ge-C (5% C) film grown on a Ge(100) substrate. 
The Ge substrate was prepared by cleaning in deionized wa¬ 
ter followed by oxide desorption in the MBE chamber at 
— 500 °C. In addition, a — 20-nm-thick Ge buffer was grown 
at 200 °C to give a sharp (2X1) reconstructed Ge surface. 

Selected samples were annealed at 750 °C for about 1 h 
in an Ar ambient to evaluate the stability of the thin films. 
Most of the thin films were characterized ex situ by trans¬ 
mission electron microscopy (TEM) and x-ray diffraction 
(XRD). The x-ray 26 scans were performed in a Siemens 
D500 diffractometer system using Cu Ka radiation. TEM 
samples were prepared by conventional mechanical polish¬ 
ing followed by ion milling at liquid nitrogen temperatures. 
The microscopy was performed in JEOL lOOcx or 4000fx 
microscopes. 

Raman spectroscopy was also performed on most 
samples. The Raman spectra were obtained with a Renishaw 
1000 Raman microscope using the 633 nm HeNe laser for 
excitation. This instrument has holographic optics to reject 
the Rayleigh scattering, a CCD array detector, and a single 
stage spectrometer. The laser power was typically 4 to 5 
mW; a 50x objective was used, giving a spot size of a few 
/urn’s; exposure times were typically 5 min, giving —10 5 
counts on the strongest lines; and the spectral resolution [full 
width at half maximum (FWHM)] was 4.7 cm -1 , as mea¬ 
sured on a series of plasma emission lines used for calibra¬ 
tion. There is sufficient transmission in our films at 633 nm 
so that all spectra showed the sharp c-Si line near 520 cm” 1 
from the substrate. This line gave a convenient reference to 
check the reproducibility of our frequency-shift measure¬ 
ments. Over a period of several days, the c-Si line from a 
dozen samples was found to occur at 521.30±0.08 cm” 1 , 
where the uncertainty is one standard deviation. The c-Ge 
pe°k measured at the same time from a Ge wafer was at 
301.70 cm” 1 . 



FIG. 1. RHEED patterns showing: (a) 3D island formation for Ge-C(5%) 
films on Si, (b) pronounced faceting of the islands (in the first few mono- 
layers) when Ge is deposited on Ge-C(5%) layer, (c) continued Ge depo¬ 
sition leading to a smooth 2D layer, and (d) tendency to form twinning 
defects when Si is deposited on Ge-C(5%). 


III. RESULTS 
A. RHEED patterns 

For Gej _ X C X thin films grown on Si, the observed evo¬ 
lution of RHEED patterns depends on the carbon content, as 
reported previously. 14 The general tendency is that at low C 
concentrations (jc< 2 to 3 at. %), the (2X1) reconstruction 
is stable and a streaky pattern indicative of a 2D layer growth 
is observed on the surface. For medium C films (x 
= 3-5 at. %), the initial 8-10 ML’s (1 ML=1.4 A) show a 
streaky (2X1) pattern, followed by evidence of transmission 
spots indicative of 3D island growth. For even higher C 
films, extra spots indicative of twinning were also observed. 

Figure l(a)-(d) shows RHEED patterns observed during 
the growth of samples A and B. Figure 1(a) is the typical 
RHEED pattern of Ge-C (5% C) films grown on Si(100). It 
is noted that the pattern shows transmission spots indicative 
of 3D growth. In addition, the spots show a chevron making 
a maximum angle with the {100} plane of —10°. Figure 1(b) 
is also a pattern from sample A, but after the first few MLs of 
pure Ge have been deposited on the Ge-C surface. Strong 
streaks indicative of facet formation on the 3D islands were 
immediately observed. The angle between the main {100} 
and the facet plane is —20°, consistent with {411} type fac¬ 
ets. With continued Ge deposition, this feature gradually 
changed to a sharp (2X1) streaky pattern, as shown in Fig. 
1(c). For sample B, however, when Si was deposited on the 
Ge-C layer at 200 °C, the RHEED pattern continued to 
show the 3D transmission spots and further developed some 
extra spots indicating the formation of twins, as shown in 
Fig. 1(d). For sample C, Ge-C thin films on Ge(100) sub¬ 
strates, strong faceting developed after the first few mono- 
layers of codeposition of Ge and C, and became more pro- 


3288 J. Appl. Phys., Vo!. 82, No. 7, 1 October 1997 


Yang, Krishnamurthy, and Weber 






1.0450 


. 1.0440 


o' 

U1-0430 \- 

i 

Ql.0420 f- 


© 1.0410 - 
cS 1.0400 - i 
^1.0390 j- 


r 1.0380 {- 


1.0370 


0.00 


V 


- 


• As-grown 

-Relaxed Ge 

|X Annealed 


2.00 4.00 6.00 8.00 

Carbon Content (at.%) 


t 

i 

10.00 


FIG. 2. The d ^ ratio of the Ge-C films to the Si substrates obtained from 
XRD, plotted against the C concentration. Notice the decrease of lattice 
parameter with increased C content for as-grown films. After anneal, there is 
some decrease in lattice parameter for C<3% alloys and some increase for 
C=5%-10% alloy films relative to the as-grown films. The dashed line 
indicates the value for relaxed bulk Ge. 


nounced upon further deposition. The facet streaks were 
measured to be from {311} type planes . 15 

B. X-ray diffraction results 

X-ray diffraction 26 scans were performed on all the 
as-grown and annealed samples. To minimize the systematic 
errors, the Ge-C(400) peak positions were normalized with 
that of the Si(400) substrate peak position. In other words, 
the ratio of ^ of Ge-C to Si was calculated from the peak 
positions. This d ratio is plotted as a function of carbon 
concentration, as shown in Fig. 2. The dashed line indicates 
the ideal ratio of the relaxed bulk Ge to that of bulk Si. It can 
be noticed that all the as-grown samples have a decreasing 
lattice spacing with increasing C concentration. The ^400 of 
pure Ge thin film shows a higher value than the relaxed bulk 
Ge, indicating that there is some epitaxial strain remaining in 
the pure Ge film. It is estimated from this d ratio that ~97% 
(of the total 4%) mismatch strain has been accommodated. 
For the Ge-C films with nominal C<2 at. %, a reduction in 
the J 400 values (towards the bulk Ge value) indicates that the 
strain is relaxed by substitutional C and/or more efficient 
generation of misfit defects. For films with even higher 
nominal C content, the ^400 value reduces below that of bulk 
Ge value and is clear evidence for the existence of substitu¬ 
tional C. The lowest d ratio measured indicates that a maxi¬ 
mum of ~ 1 at. % C has been substitutionally alloyed with 
Ge (assuming Vegard’s law). It is possible that interstitial C 
may also change the d ratio . 16 

After the anneal, it is observed that while the samples 
with lower C (<3 at. %) content show a decrease in the d 
ratio, the higher C content (>~5 at. %) films show an in¬ 
crease in the d ratio relative to the as-grown values. Possible 
reasons for this difference will be discussed later in the 
article. 

C. Transmission electron microscopy results 

Plan-view and cross-secfion transmission electron mi¬ 
croscopy have been performed on most of the thin film 


samples. Some results have been previously published . 14,15 
Typically, the cross-sectional TEM images show that the 
films are of good epitaxial quality for most of the films stud¬ 
ied. Figure 3(a) is a cross-sectional image of a pure Ge film 
on Si(100), showing a majority of the misfit defects localized 
near the interface, a few threading dislocations, and no pla¬ 
nar defects. With increasing C concentration, in addition to 
the misfit defects near the interface, the planar defect (e.g., 
stacking faults and microtwins) density is increased, as seen 
from Figs. 3(b) and 3(c). Furthermore, the growth fronts 
change from being flat to microscopically rough with in¬ 
creasing C contents, which is consistent with our RHEED 
observations. 

Plan-view TEM images show that the pure Ge and 
Ge-C films on Si substrates are relaxed with short interfacial 
misfit dislocations as evidenced by the presence of moire 
fringes and from weak-beam dark-field images as shown, for 
example, in Fig. 3(d). 

After annealing, there is a major difference between the 
low carbon content (<5 at. %) films and the higher carbon 
content films. For lower carbon content films, the misfit dis¬ 
location density seems to have increased while the threading 
dislocation and the planar defect density seem to have de¬ 
creased relative to the as-grown films as seen by comparing 
Fig. 3(a) with 4(a) and Fig. 3(b) with 4(b). For higher C 
films, the most significant change appears to have occurred 
at the interface as seen by comparing Fig. 3(c) to 4(c). A 2 to 
3 -nm-thick, bright interfacial layer has developed where the 
lattice fringes are not resolved and is likely to be amorphous. 
On the substrate side of the interfacial layer, a darker layer is 
also apparent. It is likely that the bright amorphous layer is C 
and/or Si rich (lower atomic mass leading to higher electron 
transmission), while the darker contrast on the substrate side 
is from Ge enrichment. The possible origins of the amor¬ 
phous layer will be addressed in the next section. 

Figure 5(a) is a cross-sectional TEM image of sample C, 
Ge-C thin film on Ge substrate. The image shows several 
interesting features. First, the growth front shows sharp {311} 
faceting. Second, within the thin film layer, dark lines of 
contrast making an angle consistent with {311} type planes is 
observed. This contrast is likely due to strain associated with 
C atoms that appear to have preferentially segregated to 
{311} planes. Third, unlike growth on Si substrates, there are 
no signs of misfit or planar defects. The {311} faceting and 
roughening is consistent with the RHEED patterns observed 
for this sample. The plan-view image of this sample, Fig. 
5 (b), shows no moire fringes but some strain contrast possi¬ 
bly due to C clusters. 

D. Raman spectroscopy studies 

Raman spectra of S\\ -^Ge* alloys have been studied in 
detail both experimentally 17 ” 21 and theoretically . 19,22 ” 25 
These spectra show lines near 300, 400, and 500 cm” 1 , that 
correspond, respectively, to the Ge-Ge, Ge-Si, and Si-Si 
local vibrational modes. The frequencies of these modes in 
free-standing alloys shift in systematic ways with alloy com¬ 
position according to the combined effects of mass disorder 
and microscopic strain . 23 Similar effects are seen in Sii_ x C x 
and Si-rich Si^ x -yGc x C y alloys, where the presence of a 
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FIG. 3. (a) XTEM of Ge/Si, showing mostly relaxed film with some misfit dislocations at the interface, (b) XTEM of as-grown Ge-C (C= 5%) film, showing 
planar defects in the film, (c) XTEM of as-grown Ge-C (O 10%) film with even larger planar defect density and some polycrystals developing, and (d) 
Flan-view of TEM (weak beam, dark field) of Ge-C (C=5%) film showing short misfit inf facial dislocations, similar to pure Ge on Si. 


small amount (~1%) of substitutional C in the lattice is 
manifested by the appearance of an additional line near 
600 cm -1 associated with the Si-C local mode . 7,23,26-28 In 
the Gej.jC* alloys of the present study, we see no clear 
evidence for Ge-C local mode to confirm the existence of 
substitutional C. This signal is possibly masked by a weak 
second order Ge peak that occurs near the frequency range of 
interest. There are, however, systematic C-induced changes 
in the band gap , 29 in the lattice constants and the optical 
properties 14 and similar changes in the Ge-Ge mode fre¬ 
quency, to be discussed in more detail here. We find that C 
induces a small negative shift of the Ge-Ge mode in the 
binary alloys, which we attribute to strain relaxation. As the 
C concentration increases, the Ge-Ge Raman mode develops 
a low-frequency shoulder, which we show is consistent with 
the Raman signature from amorphous Ge. In addition, an 
amorphous C peak is seen from all samples with nominal 
C>~2 to 3 at. %. 

The first question we wish to address in this section is 
the possibility that a Ge-C local mode can be seen in the 


Raman spectra. Melendez-Lira et ai 2S suggest that a feature 
they observed near 560 cm -1 in Si } _j-yGe^Cy alloys could 
be assigned to a Ge-C mode or to a Si-C mode with two C 
atoms as second neighbors. This is indeed the frequency re¬ 
gion in which one would expect to find the Ge-C mode. We 
see a similar feature in the Ge-C binary alloys, for which the 
possibility of a Si-C vibration is ruled out. However, as 
discussed below, this feature appears to coincide with a weak 
second-order Ge peak. 

Figure 6 shows three spectra obtained under identical 
conditions from a Ge(100) wafer (lower curve), from a pure 
Ge epilayer on Si(100) (middle curve), and from a nominal 
Geo. 93 Co .07 epilayer on Si(100). The Ge-wafer spectrum 
shows the strong Ge-Ge mode at 301.7 cm -1 and a number 
of weak second-order features, despite the fact that the po¬ 
larization geometry has been chosen to allow first-order and 
suppress second-order scattering. The spectra from the films 
show these same lines plus the c-Si line at 521 cm -1 and the 
second-order Si band above 900 cm -1 . The region of interest 
is the valley between the second-order Ge peak at 575 cm -1 
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FIG. 5. (a) XTEM of Ge-C(5% C) on Ge. Note the strong (311) faceting 
and C dusters decorating {311} planes, (b) Plan-view of Ge-C(5% C) on Ge 
shows strong strain contrast but no moire fringes. 



FIG. 4. (a) XTEM of Ge on Si after annealing showing significantly reduced 
threading defects, (b) XTEM of Ge-C (5% C) on Si after annealing, show¬ 
ing more misfit dislocations on the interface and less planar defects within 
the film, and (c) XTEM of Ge-C (10% C) after annealing. Note the white 
straight amorphous layer near the interface, possibly due to formation of 
a-C or silicon carbide. 


and the c -Si peak. Note that this valley is less deep in the 
C-doped film than it is in the pure-Ge film, suggesting that 
there may be some extra intensity in this region due to the 
Ge-C local mode. However, the overlap from the wings of 


the Si peak and the second-order Ge peak make it impossible 
to extract a “peak” attributable to a Ge-C local mode. Fur¬ 
thermore there is extra intensity, probably fluorescence, in 
the entire spectrum from the C-doped film. In other words, 
there is no clear signature in the Raman spectra of such a 
local mode. With a shorter wavelength pump laser, which 
would not penetrate to the substrate, o** with much thicker 
films, the interference from the Si peak’could be eliminated, 
but the second-order Ge feature will always interfere. 

The shift in the frequency of the Ge-Ge Raman mode 
produced by C doping is extremely small, typically less than 
10%-20% of the line width, and less than the frequency in¬ 
terval between points in the spectra. In order to measure this 
shift accurately, it is imperative that we use an analytical 
lineshape function that gives a good fit to the experimental 
data. The peak frequency will then correspond to the maxi¬ 
mum value of this function. We have found that simple 
Lorentzians or Gaussians do not give a good fit to the ex¬ 
perimental lineshapes, even for the pure Ge films. Moreover, 
with C added, the line becomes asymmetric. We have thus 
chosen to use Pearson-type lineshape functions to fit these 
lines. 30 For the pure Ge films, we use a symmetrical Pearson 
VII distribution: 

j(o— tfi \ 2 -ff3 

/(«) = *<> 1+4 ("^"J < 2 3 ” !) • 0) 

where a 0 is the amplitude, a x is the peak frequency, a 2 is 
the width, and a 3 is a shape parameter. When a 3 = 1, the 
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FIG. 7. Data points are polarized Raman spectra, as in Fig. 6, of the Ge-Ge 
mode, after subtracting a linear background, for three epitaxial films grown 
on Si(100). The solid lines are least-squares fits to Eq. (1) for the pure Ge 
film and to Eq. (2) for the Ge-C films. The spectra have been scaled slightly 
to have the same peak amplitude. 


FIG. 6. Polarized Raman spectra from an oriented wafer of Ge (100) (lower 
curve), a pure Ge film on Si(100) (middle curve), and a Ge-C(7%) film on 
Si(100) (upper curve). 


function is a Lorentzian, and as a 3 becomes large the shape 
tends toward a Gaussian. This lineshape is widely used as a 
simple analytical approximation to the better-known Voight 
function. 

For the C-doped films, we use an asymmetric Pearson IV 
lineshape: 


/(<o) = a 0 (l + * 2 ) “ 3 exp 



(tan 'xj + tan* 1 



(2) 


with 

cd — a\ — a 2 a 4 l2a^ 

x - -, 

a 2 

where a 0 is the amplitude, a x is the peak frequency, a 2 is a 
width, and ar 3 and a 4 are shape parameters. The asymmetry 
is determined by the sign of a 4 , a positive value gives a 
low-frequency tail, a negative value a high-frequency tail. 

Examples of our spectra along with fits to Eqs. (1) and 
(2) are shown in Fig. 7. The data points here are the raw 
spectra from which a linear background has been subtracted. 
For plotting purposes, the data for two of the plots have then 
been scaled slightly to yield the same peak intensity as the 
third. The solid lines in this figure are least-squares fits to the 
Pearson functions given above. Note that there is a small 
negative shift of the peak frequency with increasing C and 
that there is a substantial tail on the low-frequency side of 


the 5% C curve. A summary of all the frequency-shift data is 
given in Fig. 8. The error bars in this figure are 2.5 times the 
standard deviation of c-Si peak frequency measured at the 
same time from the same spectra. We consider these error 
bars to be conservative estimates of the experimental uncer- 



FIG. 8. Plot of the Ge-Ge frequency, determined by fits as in Fig. 7, as a 
function of C concentration in the films. The lower left point is for a Ge 
wafer. The dashed line is a linear least-squares fit to show the trend. 
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FIG. 9. Raman spectra after linear background substraction of amorphous 
(a) and crystalline (b) Ge-C films. The data points and the resultant fits 
(solid curves) are shifted up by one vertical tick mark for clarity. The de¬ 
composition of the spectrum in (a) into two peaks, one for a-Gt the other 
for c-Ge shown below the data, is done with a full least-squares minimiza¬ 
tion. The a-Ge lineshape determined in (a) is then constrained to be the 
same in (b). 



FIG. 10. Unpolarized Raman spectra from Ge-C films showing the broad 
band from amorphous C at the highest concentrations. The spectra have 
been systematically shifted by a few hundred counts in order to make a 
nearly evenly spaced stacked plot. 


tainty, since the Ge peak is actually bigger than the c-Si 
peak. Although there is considerable scatter in the points, 
the trend toward lower frequency with increasing C is 
convincing. 

There are three effects contributing to the shifts shown 
in Fig. 8: mass disorder, microscopic strain, and epitaxial 
strain due to the lattice mismatch between film and substrate. 
The last of these can be calculated from measurements of the 
film lattice constants and the known deformation parameters 
for the phonon. Our x-ray data indicate that the pure Ge films 
have positive strain in the direction normal to the surface, 
£ zz = 2.4X 10“\ and a corresponding negative in-plane 
strain, e xx = — (C u /2Ci 2 )€ zz = — 3.2X 10“ 3 , where the C (; 
are the Ge elastic constants. This biaxial strain will split the 
triply degenerate zone-centered Raman phonon into a singlet 
and a doublet. Only the singlet will be observed in our scat¬ 
tering geometry, and its dependence on the strain is given by 
Refs. 23 and 31: 


(O— co 0 


p — q 

l-r(2e„+ e zz ) - e zz ) 


(3) 


where a> 0 is the frequency in an unstrained lattice, y is the 
Griineisen parameter, and (p — q)/ 3o)q is a coefficient de¬ 
scribing the strain dependence of the internal force field. Us¬ 
ing the values y=0.89 and (p — q)/3^0 = 0.15, measured by 
Cerdeira et al ., 31 we calculate a shift of +1.3 cm -1 for the 


epitaxial film. This value agrees well with our measured 
value of 1.2±0.2 cm” 1 from the data in Fig. 8. 

As the C concentration increases, the Ge lattice param¬ 
eter normal to the surface decreases, which has the tendency 
to relax the strain (island formation is also relaxing the 
strain, but it is difficult to separate the two effects). The 
strain relaxation shifts the peak towards the frequency for 
bulk Ge. Based on the x-ray diffraction data, we notice that 
the strain is already relaxed (i.e., bulk Ge lattice parameter is 
reached) for films with more than ~ 2 at. % nominal C. 
However, the Raman data still show a positive shift com¬ 
pared to bulk Ge. It is therefore likely that this shift is a net 
result of the influence of C-induced microscopic strain and 
mass disorder. 

A third interesting feature observed in our Raman data is 
that at the highest C concentrations a broad shoulder near 
280 cm” 1 is observed that cannot be fit with a single Pearson 
IV lineshape. We accidentally discovered an extreme ex¬ 
ample of this shoulder in an attempt to grow a 6 at. % C film 
on a substrate from which the oxide was not completely 
removed. The result was a nearly amorphous film, as indi¬ 
cated by x-ray analysis. The Raman spectrum from this film 
is shown in Fig. 9(a). This spectrum is composed of two 
peaks, a broad highly asymmetric band centered near 
280 cm” 1 and a narrow peak at 301.6 cm” 1 . The broad band 
is the signature for amorphous Ge, 32,33 and the narrow peak 
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FIG. II. Amorphous C Raman peak intensity as a function of C concentra¬ 
tion as described in the text. The straight line is a linear least-squares fit. 

is from Ge crystallites embedded in the mostly amorphous 
film. The frequency of the narrow peak matches well with 
that for fully relaxed crystalline Ge. 

The spectrum in Fig. 9(b) is from an epitaxial Ge 09 Co.i 
film. It is fit by using the same lineshape for the broad band 
from Fig. 9(a), in which only the amplitude is varied, along 
with the usual Pearson IV lineshape for the narrow peak. The 
fit is convincing. The broad shoulder is completely ac¬ 
counted for by allowing for a small contribution from the 
a-Gc lineshape. We conclude that the high C films have an 
increasing fraction of amorphous Ge. 

Raman measurements were also made in the 
800-1800 cm -1 range to explore the possibility of C-related 
vibrational modes. Figure 10 shows a series of spectra for 
various Ge-C alloys. The signature of amorphous C, a broad 
band peaking near 1400 cm -1 is clearly seen lor the higher C 
films. 34 The 10 at. % C film was fitted to gtl the lineshape, 
which was then constrained to obtain the peak strength for 
the rest of the plots, as shown in Fig. 11. A straight line fit to 
the data gives an intercept of 2.4% below which no a-C peak 
is present. 

Annealing produces a major change in the C-C region 
of the spectrum, as shown in Fig. 12. The annealed films 
have a double-peak structure characteristic of graphitic car¬ 
bon. There is a fairly broad line near 1330 cm" 1 , corre¬ 
sponding to the disorder-induced zone-boundary phonon 
from graphite and a narrower line at 1598 cm -1 , correspond¬ 
ing to the zone-center E 2g mode. This material is also called 
glassy carbon or microcrystalline graphite. 35,36 The graphitic 
C bands are much stronger than the a-C band in the unan¬ 
nealed samples, but that is probably due to the resonance 
effect owing to a strong optical absorption in graphite. 

IV. DISCUSSION 

A. Incorporation of C in Ge^C* films 

Our experimental results ane consistent with the follow¬ 
ing model. When the nominal C concentration is less than 



Raman shift (cm" 1 ) 


FIG. 12. Effect of annealing on a nominal Ge-C (10%) film. Note the 
change in the a-C peak to a graphitic C double-peaked structure. 

~2 to 3 at. %, the deposited C is predominantly dispersed 
either substitutionally and/or interstitially. This results in the 
same growth mode as pure Ge, aside from the fact that the 
films have a lower lattice parameter due to substitutional 
incorporation of C. When the increase in the C flux leads to 
an increase in the nominal C>~3 at. %, not only does more 
C go into Ge substitutionally (as indicated by the x-ray data), 
but also the growth mode changes from 2D to 3D islanding/ 
facet formation after initial few MLs of Ge-C codeposition. 
While there may be a marginal increase in the substitutional 
fraction due to an increase in the C growth rate, 8 it is sug¬ 
gested, based on correlation of the a-C observed by Raman 
with the RHEED and TEM results, that the excess a-C is 
responsible for the change in growth mode from 2D layer to 
3D island growth (discussed in Sec. IV B). The formation of 
planar defects in the Gej_ x C^ films with nominal 03 at. % 
is likely to be a consequence of island coalescence in the 
presence of strain (In the absence of significant strain, e.g., 
Ge-C on Ge, planar defects do not seem to form). Conse¬ 
quently, the tendency for island formation increases with in¬ 
creased nominal C concentration leading to an increased 
density of planar defects. 

In the particular case of the G^i- x C x (x—5 at. %) films 
grown on Ge substrates, the epitaxial mismatch strain is un¬ 
likely to be very high, based on the fact that no extra Ge-C 
peak can be detected in the x-ray 26 scans. This is also 
consistent with our observation that no misfit induced defects 
are found in TEM studies. In this case, the excess C not only 
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changes the growth mode from layer to 3D islands, but also 
the excess C is forced to precipitate/segregate along {311} 
planes, as seen in Fig. 5(a). One likely reason as to why the 
{311} planes are specifically decorated in Ge^C* films 
grown on Ge may be related to a lower surface energy of 
Ge-C {311} compared to Ge-C {100} planes, as calculated 
for Si surfaces. 37,38 In the case of Ge^^C* films on Si, the 
tendency to facet along {311} is less pronounced possibly 
due to kinetic constraints (growth of Ge 0 . 95 C 00 5 on Si at 
550 °C leads to strong {311} faceting). The lack of decora¬ 
tion of {311} facets in Ge^C^/Si may be due to the fact 
that some of the C may be incorporated at the high density of 
defect (planar and misfit) sites and from epitaxial strain ef¬ 
fects. 

B. The driving force for 2D to 3D transition 

It is well known that pure Ge grown on Si follows the 
layer plus island or the Stranski-Krastanow growth mode. 
This is because island formation relieves the strain energy at 
the expense of increased surface area. However, by maintain¬ 
ing a sufficiently low substrate temperature (e.g., 200 °C), 
the reduced surface diffusion suppresses island formation, as 
confirmed by our RHEED and TEM observations. In our 
studies of Ge-C alloys, the presence of 03 at. % leads to 
3D island formation even at this low substrate temperature 
(200 °C). The possible reasons 14 for this behavior could be 
(i) increased strain-induced surface diffusion in the presence 
of C, 39 (ii) changes in surface energetics/wetting behavior 
due to C, or (iii) heterogeneous nucleation at C clusters. 40 

We rule out nucleation at C clusters as a possibility be¬ 
cause Ge growth on a rough Ge-C surface rapidly becomes 
planar at this growth temperature. While increased surface 
diffusion is a possibility, the islanding tendency is not due to 
strain alone as evident from the fact that Ge-C grown on Ge 
tends to grow 3D. The most likely reason appears to be a 
change in surface energies of the Ge-C film relative to Si 
and Ge, leading to either changes in the wetting behavior or 
an increased tendency to lower their surface energy by form¬ 
ing specific crystallographic facets, e.g., {311}. 

The wetting behavior may be explained on the basis of 
the following qualitative argument. It is proposed that the 
surface energy of Ge-C is larger than that of Ge and is a 
linear function of C concentration (It is assumed that the 
interfacial energies involved are all the same). The surface 
energy of pure Si is known to be larger than that of pure Ge 
and explains why Ge wets Si. On that basis, it can be under¬ 
stood that a Ge-C film with low C concentration would have 
a surface energy lower than that of Si and is likely to grow in 
a layer mode. Ge-C films with higher C concentration and 
whose surface energy is larger than that of Si may grow in an 
island mode. In our experiments, this cross-over appears at a 
C concentration of about 3 at. %. This C rich surface subse¬ 
quently lowers its energy by faceting along {311} planes. 
This model explains the observed morphology in the experi¬ 
ments involving the growth of Ge-C on Si, Ge-C on Ge, 
and Ge on Ge-C films at 200 °C. The only apparent excep¬ 
tion is the growth of Si on Ge^Q) 05 where 3D growth was 
observed after 20 nm of Si deposition, even though our 
model would predict 2D growth. It is however conceivable 


that by forming facets, the Ge-C may have reduced its sur¬ 
face energy to a value below that of Si. An alternative ex¬ 
planation is the lack of sufficient kinetics for Si to smooth 
out the roughness on a Ge-C (100) surface. 

C. Effect of annealing 

Annealing at 750 °C causes measurable changes in the 
microstructure and composition of the Ge-C alloys. From 
the x-ray data (Fig. 2) for the annealed samples, it can be 
noted that for pure Ge and low C(<3 at. %) concentration 
alloys, annealing leads to a marginal reduction in the lattice 
parameter below that of bulk Ge (i.e., the dashed line). One 
cause for the decrease in lattice parameter could be that the 
amount of substitutional C may have increased. In addition, 
Si interdiffusion into the films during the anneal may also 
contribute to this effect. The presence of ^2-4 at. % Si has 
been measured in our annealed samples by Raman 
spectroscopy 41 The postannealed samples also appear to 
have a higher density of interfacial misfit dislocations as seen 
by comparison of Figs. 3 and 4, while the density of planar 
defects and threading defects seems to have decreased sig¬ 
nificantly. All of these effects can contribute to the decrease 
of the lattice parameter relative to the as-grown films. 

For higher C films (>5 at. %), the x-ray data indicate 
that the lattice parameter increases after the anneal. A likely 
explanation is that a significant amount of C has been ejected 
out of the substitutional sites, and that this effect outweighs 
the reduction in lattice parameter due to interdiffused Si. A 
second microstructural feature of significance is the amor¬ 
phous layer formed at the Ge-C/Si interface. This layer ap¬ 
pears only after the anneal and seems to be rich in low 
atomic number elements (e.g., C) as deduced from the TEM 
contrast. Since the diffusivity of Si and Ge is likely to be 
larger than C, and a driving force exists to form the stable 
compound SiC, a C-rich Si-C amorphous layer may have 
formed a' the interface. Further work is needed to character¬ 
ize the composition and understand the formation of this 
layer. 

The incorporation of C and the growth mechanism of 
Ge-C alloys on Si at low temperature can be summarized as 
follows. We start with the assumption that there is a limited 
“substitutional solubility” and “interstitial solubility” of C 
in Ge in the thin film state. This solubility depends on the 
growth conditions (e.g., temperature, flux, and the strain). 
Since most of our films were on Si substrates and grown at 
200 °C, the growth temperature and the epitaxial strain are 
maintained nearly constant. Under these conditions, the equi¬ 
librium solubilities are almost constant (depending only on 
the C flux). At low nominal C content, the deposited C is 
distributed both substitutionally and interstitially and the 
growth front is planar. Beyond the solubility limit, the excess 
C exists in the form of a~C on the surface (i.e., at the growth 
front) and alters the growth mode from 2D layer to 3D 
growth. This transition occurs about 3 at. % in our studies. 
For C contents larger than 5 at. %, the coalescence of the 
3D islands (in the presence of strain) leads to the formation 
of a high concentration of planar defects (e.g., stacking faults 
and microtwins) as seen from the cross-sectional TEM im- 
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ages. For even higher C contents, the presence of a high 
concentration of a-C degrades the epitaxial quality of the 
film. 
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Low-temperature (—200 °C) molecular beam epitaxy of Ge-rich Gei_ x _ySi v C x alloys grown on 
Si(100) have been investigated by in situ reflection high-energy electron diffraction, ex situ x-ray 
diffraction, transmission electron microscopy, Raman scattering, and ellipsometry. The Si contents 
were either —20 or —40 at % and the C concentrations were nominally varied from zero up to —8 
at %. Selected samples were annealed in an Ar ambient at 750 °C to evaluate the stability of the thin 
films. With increasing C concentration, the epitaxial growth mode changes from two-dimensional 
(2D) layer growth to 3D island growth. Under the growth conditions studied, the GeSiC films have 
a tendency to form planar defects, whose density increases with increasing C and Si concentrations. 
The x-ray diffraction data show that the lattice parameter decreases with increasing C concentration. 
It is estimated that a maximum of —2-3 at % C is substitutionally incorporated into these films. 
Raman spectra of the alloy films show that the effects of C on the strong Ge-Ge and Ge-Si local 
modes are far less than the effects due to Si. We are unable to observe any systematic change in the 
Ge-Ge mode, whereas the Ge-Si mode appears to shift to lower frequency with the small addition 
of C. Ge 1 ^ x .. } ,Si > ,C x films formed by annealing Gej_ x C x films on Si are also discussed. 
Spectroscopic ellipsometry determinations of the film’s optical constants show that the primary 
effect of C is to reduce the strength of the E x critical point feature. © 1998 American Institute of 
Physics. [S0021 -8979(98)02316-0] 


C levels up to —2-3 at %, 6,10 with little emphasis on the 
Ge-rich ternary system, except for Ref. 10. For a review of 
the growth and properties of SiGeC alloys, see Ref. 13. 

We have recently reported 14,15 on the microstructural 
evolution and optical properties of the Gt\- x C x alloys 
(x<0.1) grown at —200 °C. In this article we present the 
effect of adding Si to Gej _ X C X alloys grown on Si(100) sub¬ 
strates at —200 °C. The reasons for using such a low tem¬ 
perature (reduced kinetics) are to suppress phase separation 
and to prevent SiC precipitation. In addition, it has also been 
observed in the Si-C system that a reduction in growth tem¬ 
perature leads to an increase in the substitutional-to- 
interstitial carbon ratio. 16 Our selection of 200 °C is based on 
the fact that it is an estimated lower limit for good epitaxial 
growth. The films were characterized in situ by reflection 
high-energy electron diffraction (RHEED) and ex situ by 
transmission electron microscopy (TEM) and x-ray diffrac¬ 
tion (XRD). Optical characterization was performed using 
Raman spectroscopy and spectroscopic ellipsometry (SE). 

II. EXPERIMENTS 

The Gej _ x _ y Si x C y films were grown in a Riber 32 MBE 
system with a base pressure of better than 5X 10” 10 Torr. 
Si(100) substrates were chemically cleaned and oxidized by 
boiling in a H 2 02 :H 2 S 04 mixture. After the sample had been 
degassed in the MBE at -200 °C overnight, the oxide was 
desorbed by heating to just below 900 °C for about 15 min 


I. INTRODUCTION 

SiGe alloys have been extensively studied in order to 
exploit the promise of band-gap engineering on silicon. 1 In 
the past few years, alloys of group IV elements including C, 
Si, Ge, and Sn have been actively investigated for use in 
heterojunction devices compatible with current Si integrated 
circuit technology. 2,3 There has been a particularly strong 
interest in SiGe alloys containing carbon. The main reason is 
that (assuming Vegard’s law) an approximately 9:1 Ge to C 
ratio is expected to lattice match SiGeC alloys to Si, circum¬ 
venting critical thickness limitations inherent in SiGe tech¬ 
nology. Furthermore, C in the SiGe lattice makes it possible 
to use both alloy concentration and strain as a variable for 
band-gap engineering. However, several obstacles must be 
overcome. First, C is relatively insoluble in both Ge and Si, 
unlike the completely miscible SiGe system. Second, the for¬ 
mation of various forms of SiC is thermodynamically fa¬ 
vored at higher growth temperatures. 

Despite these problems, low concentrations of C in Si- 
Ge-C have been successfully incorporated using a variety of 
techniques including chemical vapor deposition (CVD) , 4 ’ 8 
molecular beam epitaxy (MBE), 9 ” 11 and solid phase 
epitaxy. 12 Those thin film alloys were typically grown at 
temperatures between 400 and 650 °C, with a Ge content of 
20 at %, 3 30 at %, 5 40 at %, 6 or even less, with substitutional 
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and then the sample temperature was quickly dropped to 
—200 °C. Sharp (2X1) diffraction spots indicative of a clean 
surface were typically observed in RHEED. In all cases, in 
order to produce a better quality surface, a couple of mono- 
layers of Ge were deposited on the sample surface at 200 °C 
soon after the above oxide desorption, followed by re¬ 
heating to —800 °C for complete oxide desorption. After the 
above Ge-assisted oxide desorption, the RHEED showed a 
sharper (2X1) reconstruction pattern indicating a better qual¬ 
ity surface. Germanium was deposited from a previously 
calibrated pyrolyP', boron nitride Knudsen source at —0.5 or 
1 nm/ntin. The C was deposited from a rod-fed electron- 
beam source, and the flux was varied to obtain various C 
concentrations. Si was also deposited from a rod-fed 
electron-beam source whose flux was kept constant at —0.3 
nm/min. Rutherford backscattering spectrometry (RBS) (in¬ 
cluding C resonance) was used to estimate the C and Si 
concentrations. A typical experiment consisted of the follow¬ 
ing steps. After the sample had stabilized at —200 °C, the Ge 
shutter was first opened followed by the Si and C shutters 
after a 10-30 s delay. The nominal thickness for most of 
films varied from 30 to 70 nm. Two groups of ternary alloys 
have been studied, namely (i) type A: high Ge content alloys, 
with —75% Ge, the balance being Si about 20% and C vary¬ 
ing; (ii) type B: low Ge content alloys with —55% Ge, the 
balance being Si about —40% and C varying. 

The Raman measurements were done using a Renishaw 
1000 Raman microscope with 633 nm excitation. The experi¬ 
mental apparatus and the line shape fitting procedure are de¬ 
scribed in more detail elsewhere. 15 The spectroscopic ellip- 
sometry data were recorded with a variable-angle, automated 
ellipsometer made by J. A. Woollam Co. Optical constants 
of the films were obtained by fitting to an expansion of e as 
a sum of Lorentz oscillators representing the various critical 
points in the opucal response of Si-Ge alloys. The position, 
strength, and width *of all oscillators in the 1-5 eV range 
were optimized as were the film and native oxide thick¬ 
nesses. 


III. RESULTS 

A. Reflection high-energy electron diffraction 
(RHEED) studies 

The RHEED patterns of Ge 0 . 8 Sio. 2 , Ge^S^ and 
Geo. 7 gSio. 2 Co .02 all °y fi* ms were streaky (2X1) patterns 
that remained almost unchanged throughout the growth. On 
addition of more C (>3-4 at %) the layers show initial layer 
growth followed by 3D island formation. For even higher C 
contents (>6 at %) the patterns changed quickly from (2X1) 
streaks to 3D transmission spots, eventually forming extra 
twinning-related spots at the end of the growth. The C- 
dependent changes in RHEED patterns associated with the 
low-temperature growth of Ge-rich GeSiC alloys are similar 
to those observed for low-temperature Ge-C alloys grown 
on Si. 14 ' 15 




FIG. 1. XTEM images of (a) GeojSioj film, showing smooth growth front 
and some misfit defects at the interface, and (b) Ge^Sio^QxtM film* show¬ 
ing rough growth front and planar defects in the film. 


B. Transmission electron microscopy (TEM) 
characterization 

The microstructural features of the high Ge and low Ge 
content GeSiC alloys are characterized by short interfacial 
misfit dislocations as observed in plan-view TEM. With in¬ 
creasing C content, an increasing density of planar defects is 
observed in cross-sectional TEM (XTEM) images. Figure 1 
shows a comparison between the XTEM images of high Ge 
content GeSi alloys and GeSiC alloys. The Ge 0 gSi 0 .2 alloy 
shown in Fig. 1(a) has some misfit dislocations at the inter¬ 
face and almost no planar defects within the thin film. With 
the addition of C, a significant increase in the number of 
planar defects (mainly twins and/or stacking faults) are ob¬ 
served, as seen in Fig. 1 (b) for a Ge 0 76 Si<) 2 C 0 04 alloy film. In 
addition, the growth front of the GeSiC film is seen to be 
very rough [Fig. 1(b)] as compared to the equivalent GeSi 
alloy [Fig. 1(a)], consistent with the RHEED observations. 
Examples of XTEM images of low Ge content films are 
shown in Fig. 2. Figure 2(a) is an image from a Ge 06 Si 04 
alloy, while Figs. 2(b) and 2(c) are the images for 
Ge 0 . 59 Sio 34 C 0.7 alloys with two different film thicknesses (30 
and 70 nm). The growth front appears rough for the GeSiC 
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FIG. 2. XTEM images of (a) ~30 nm Geo.gSio^ film. Mo wing smooth 
growth front, (b) —30 nm Geoj^io^Qxn film, and (c) —70 nm 
Geo. 59 Sio. 34 Co .07 film, (b) and (c) show rough growth front and an increased 
density of planar defects in the films compared with (a). 


film in comparison to the GeSi film. However, the morphol¬ 
ogy of the thicker GeSiC film is not significantly different 
compared to the thinner film. 

Statistical analyses were performed on the various- 
XTEM images to obtain a rough estimate of the planar defect 
density as a function of C content, and are plotted in Fig. 3. 
The basic tendency is that the planar defect density increases 
with C content both for the Ge-C alloys and for the two 
groups of GeSiC alloys. The addition of Si appears to en¬ 
hance the tendency for planar defect formation. This is evi¬ 
dent from a comparison between Ge-Si-C alloy with —4-7 
at % C whose defect density is similar to Ge-C alloys with 
10 at % C . 15 This implies that the addition of 20 at % or 40 
at % Si does not improve the film quality under the condi¬ 
tions studied. 

C. X-ray diffraction (XRD) studies 

X-ray 20 scans were performed in the (400) range for all 
the samples. The ^400 ratio of the thin film alloy peak to that 



FIG. 3. Plot of planar defect density vs C concentration for Ge-C, type A 
(high Ge content) GeSiC and type B (low Ge content) GeSiC films. The 
rough estimation of the planar defect density in the films was obtained from 
XTEM images and is in the units of cm -2 . 


of the Si substrate is plotted as a function of nominal C 
content in Fig. 4. The XRD data indicate that with increasing 
C content, the ^400 spacing of the alloy films decreases, 
which suggests that C relaxes some of the compressive strain 
either by defect formation and/or by substitutional incorpo¬ 
ration. If the equivalent GeSi thin film alloy is completely 
relaxed (assuming that the high temperature annealed sample 
is almost fully relaxed), the likely reason for GeSiC alloy 
thin film having a lower lattice parameter than GeSi is by C 
substitutional incorporation. The slope of the curve (which 
indicates the rate of substitutional C incorporation) appears 
to be similar for the high Ge content and the low Ge content 
alloys. The highest substitutional C incorporation is —2.5 
at % C (based on Vegard’s law) and corresponds to the 
thicker (70 nm) Geo. 59 Sio. 34 Co .07 alloy film * An interesting 
point is that for the same composition, the thicker film ap¬ 
pears to have a lower lattice parameter than the thinner film. 
Possible reasons will be addressed later. 
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FIG. 4. The ratio of the (400) spacing of the film to that of the Si substrate 
(4fii m ld s i) plotted as a function of nominal C concentration for both type A 
(high Ge content GeSiC) and type B (low Ge content GeSiC) films in the 
as-grown and post-annealed conditions. The single solid box and solid dia¬ 
mond data points (showing the lowest value of d^yo) are for the thicker low 
Ge content GeSiC films in the as-grown and post-annealed conditions. 
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FIG. 5. Survey spectrum of a 70 nm Geo^Si^ epitaxial film on Si(100) 
obtained with 3 mW of 633 nm laser and 300 s exposure time. 


D. Effect of annealing 

Most samples were annealed ex situ in Ar atmosphere at 
750 °C for 1 h. The XRD data (Fig. 4) show a small increase 
in the spacings compared with the corresponding 
samples before the anneal. The change is however minimal 
for the GeSi samples (especially the Geo.gSio. 2 ) suggesting 
perhaps that the as-grown sample is almost completely strain 
relaxed. It is also noted that the samples with the same com¬ 
position (Ge 0 5 9 Si{) 3 4 C() 07 ) but different thicknesses did not 
approach the same lattice parameter after the anneal. This 
result indicates that more C is incorporated in the thicker 
film than in the thinner film. 

Another microstructural difference observed after an¬ 
nealing is that the XTEM images (not shown here) show that 
the misfit defect density at the interface has increased and the 
planar defect density within the film has decreased. This is 
true for both groups of GeSiC alloys. The microstructural 
effects of annealing Ge-C films on Si have already been 
published before . 15 The annealing causes Si to diffuse into 
the Ge-C layer as discussed in the following Raman spec¬ 
troscopy section. 


E. Raman spectroscopy 

Figure 5 shows a typical unpolarized survey Raman 
spectrum of a 70 nm Ge^ g Si <).2 film. The main features in this 
spectrum have been discussed at length in previous 
studies . 17 " 25 The sharp peaks near 300, 400, and 500 cm ” 1 
arise from the Ge-Ge, Ge-Si, and c-Si first-order modes, 
respectively. The remaining structure arises primarily from 
second-order scattering, as indicated in the figure. The film is 
sufficiently transparent at 633 nm that even the weak second- 
order features from the substrate are observable. Note that 
the second-order scattering from the Ge-Ge and Ge-Si 
modes of the film can be easily seen as can the first-order 
structure near 430 and 460 cm ” 1 associated with differing 
random environments around Ge-Si and Si-Si pairs, 
respectively . 19,22 



FIG. 6. Frequencies of the Ge-Ge and Ge-Si local modes as a function of 
the Si percentage in Ge-Si-C films on Si(100). 

When C is added to a GeSi alloy we would expect sys¬ 
tematic shifts of the primary Ge-Ge and Ge-Si modes in the 
film similar to those seen in Si-rich alloys by Melendez-Lira 
et al. 26 and in Ge^ x C x alloys by Yang et al} 5 Because of 
the low concentration of Si and interference from second- 
order substrate scattering, we are unable to see the Si-C 
local mode, which has been reported at 605 cm ” 1 in 
Sij - x C x s~ Sy21 The corresponding Ge-C local mode has not 
been n^vioiisly seen, but its frequency can be estimated 
from the 605 cm ” 1 value using the following argument. The 
larger reduced mass of the Ge-C pair compared with Si-C 
will shift the frequency down by ~9.7%, and the larger Ge 
lattice constant will lead to an additional, smaller downward 
shift. Thus, the Ge-C local mode in Ge-rich alloys is ex¬ 
pected to be in the 520-540 cm ” 1 region and will likely be 
obscured in the present study by the strong c-Si peak of the 
substrate. We have recently seen this mode at 529.8 cm ” 1 in 
Gej_ x C x films grown on Ge(100) substrates, for which there 
is no interference 28 

Figures 6 (a) and 6 (b) summarize the results on the 
Ge-Ge and Ge-Si mode frequencies for a large number of 
samples plotted as a function of the Si content. The spectra 
used for these data were done in the z(xy)z” scattering ge¬ 
ometry, in order to reduce the second-order structure. The 
line frequencies in each case were determined by fitting with 
Pearson-type line shape functions, as described previously . 15 
For comparison the trends observed for these frequencies in 
bulk, unstrained, polycrystalline alloys by Renucci et al. 11 
are shown as dashed lines. 

For the Ge-Ge frequency, shown in Fig. 6 (a), there ap¬ 
pear to be no obvious C effects. However, there is sufficient 
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scatter in the data that a small shift, of order 1-2 cm” 1 as 
seen in Gci^ x C x films, 15 cannot be ruled out. On the other 
hand, large positive shifts of 6-7 cm” 1 seen in Si-rich alloys 
with fully strained, CVD-grown films can be ruled out. 26 As 
shown in Fig. 6(b) for the samples with nominal Si fractions 
of 20%-40%, C produces a negative shift of the Ge-Si fre¬ 
quency of a few cm -1 . This shift is opposite in sign and 
smaller in magnitude than the shifts seen for this mode in the 
CVD-grown film. 26 These differences are likely caused by 
differences in the strain as well as in the Ge fraction. 

We have also included in Fig. 6 the results for Ge^C* 
films into which Si has diffused during the annealing 
process. 15 The incorporation of Si in these films is demon¬ 
strated by the appearance of the Ge-Si local mode. The Si 
concentration has been estimated from the ratio of the inte¬ 
grated intensities of the Ge-Ge to Ge-Si lines. This ratio is 
expected to have a functional form A( 1 —x)!2x for a random 
alloy. We first fit this form with data from the films with 
known Si content, obtaining the value A = 2.19. Then the 
same value of A was used with the experimental intensity 
ratio to determine x for the films into which an unknown 
amount of Si had diffused. The resulting x values were 2%- 
6%. These are only rough estimates of Jt, since the Si con¬ 
centration is not uniform. However, the general trend for 
these points as far as the dependence on Si % is the same as 
for the higher Si-content films. 
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FIG. 7. Optical constants of GeSiC films on Si(100) obtained from spectro¬ 
scopic ellipsometry results. 


F. Spectroscopic ellipsometry 

Spectroscopic ellipsometry measurements were made on 
the higher Si-content films for at least three angles of inci¬ 
dence in the range 60°-80° and covering the 1-5 eV energy 
range. Results are shown in Fig. 7 for films of nominal 20% 
Si content, with and without C; similar results are shown in 
Fig. 8 for films with —40% Si. The linear increase in energy 
of the E x critical point (featured in 2-2.5 eV range) with 
increasing Si content is observed with all of our films. The 
slight energy shifts in the E\ structure in Figs. 7 and 8 can be 
attributed to differences in Si content. The primary effect of 
C is a weakening of the structure near the E x point. This 
effect is observed in all films at all C concentrations and is 
consistent with our earlier observations on Gei-jC* 
samples. 14 

IV. DISCUSSION 

A. Driving force for 3D growth 

The XTEM images of GeSi alloys [Figs. 1(a) and 2(a)] 
showed that the growth fronts were planar, while the alloys 
with C show 3D growth morphology. The smooth morphol¬ 
ogy of GeSi alloys can be understood based on the argument 
that the only driving force for 3D island formation is strain 
energy relaxation, which is suppressed by a low growth tem¬ 
perature (—200 °C). This is consistent with our previous 
study of pure Ge on Si. 14,15 However, when C is present 
beyond a certain amount, the growth mode changes from 2D 
to 3D for both Ge-C and GeSiC thin films. As discussed in 
detail previously, 14,15 it is most likely that the presence of C 
(both in crystalline form and in the form of amorphous C) 


changes the surface energy such that Ge-C and GeSiC films 
find it unfavorable to wet the surface, leading to 3D island 
formation. From our RHEED observations, it is noted that 
the larger the fraction of nominal C in the flux, the sooner 3D 
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FIG. 8. Optical constants of GeSiC films on Si(100) obtained from spectro¬ 
scopic ellipsometry results. The difference in the Ej position of the two 
samples is primarily due to their different Si content. 





2016 J. Appl. Phys., Vol. 84, No. 4, 15 August 1998 


Yang, Krishnamurthy, and Weber 


islands tend to form. So it is conceivable that a larger driving 
force for 3D island formation (proportional to the amount of 
C) exists for SiGeC alloys such that island formation is ob¬ 
served even at such low substrate temperatures. 

B. Formation of planar defects 

From our experimental observation (Fig. 3), it is clear 
that under our growth conditions, both Ge-C and GeSiC 
alloy films form planar defects unlike pure Ge films on Si. 
Furthermore, the planar defect density increases with in¬ 
creasing C and Si concentrations and is more pronounced 
when both Si and C are present. The SiGeC alloys contain a 
higher planar defect density compared with equivalent Ge-C 
and GeSi films. Several possible origins of planar defect (es¬ 
pecially stacking fault) formation have been suggested pre¬ 
viously, including: the accommodation of misfit between 
coalescing islands, 29 nucleation of faults at heterogeneous 
SiC nanocrystallites, 13 or the formation of partial disloca¬ 
tions that move towards the surface, leaving stacking faults 
behind. 30 Our results suggest that the planar defect (stacking 
faults and twins) formation is closely related to the surface 
3D morphology and therefore most likely to have formed 
from coalescence of islands. Furthermore, a clear correlation 
exists between increasing C content, 3D island density, and 
planar defect density. 

The presence of Si in Ge-C amplifies the 3D growth 
morphology. We believe that the lower mobility of Si on the 
surface enhances 3D island formation and therefore increases 
the planar defect density. Such an idea is also consistent with 
our observation that high Si content GeSi films show more 
planar defects than low Si content films [compare Figs. 1(a) 
and 2(a)]. It is interesting to observe some similarities be¬ 
tween the formation of planar defects in SiGeC films (espe¬ 
cially the effect of Si) and low-temperature Si homoepitaxial 
growth (thickness-limited epitaxy). 31,32 It has been reported 
that surface roughening leads to the formation of facets and 
twins and eventual breakdown in crystalline epitaxy. While 
the origin of 3D growth is still unclear, effects such as step- 
edge barriers have been invoked at low growth 
temperatures. 31,32 We believe that such effects may also be 
contributing to the 3D growth in our films. 

C. Incorporation of substitutional carbon 

From the x-ray data (shown in Fig. 4), it can be noted 
that the spacings of the Ge-Si-C alloys have a ten¬ 
dency to decrease with increasing C content. The decrease in 
the lattice parameter to a value lower than the almost fully 
relaxed GeSi (assuming the high temperature anneal results 
in near-complete relaxation of the GeSi alloy) suggests that 
C reduces the lattice parameter in the films by substitutional 
incorporation, which is consistent with other results, e.g.,. 4,12 
The maximum amount of substitutional C is estimated to be 
—2-3 at % based on the lowest value given in Fig. 4. This 
value is close to the data given by Refs. 6 and 10. The likely 
reason as to why the post-annealed samples have a higher 
lattice parameter than the as-grown samples, is that some 
substitutional C may have been ejected out from the films 
during high-temperature anneal. 


D. Effect of film thickness 

The interesting behavior of the samples with the same 
composition, but two different thicknesses (where the thicker 
film has a lower lattice parameter in both as-grown and an¬ 
nealed conditions) needs to be addressed. While the exact 
reason for this behavior is still unclear and further work is in 
progress, we suggest the following scenario. One possibility 
is that the thicker film is likely to have relaxed more strain 
than the thinner film. Consequently, the difference in the 
lattice parameter could be simply due to the larger residual 
compressive strain in the thinner film. Alternatively, it is 
possible that more substitutional C is incorporated during 
growth in the thicker film because of the lower strain in the 
film. We believe the latter effect is more likely because the 
lattice parameter of the post-annealed thicker film is still 
lower than the thinner film. We speculate that the C incor¬ 
poration in such films may be a function of film thickness, 
which in turn is related to strain relaxation. 

E. Effect of annealing 

The annealing experiments at 750 °C show that the 
GeSiC alloys seem to have altered the microstructure, e.g., 
decreasing the planar defects density, leading to an improve¬ 
ment in the film quality. On the other hand, some substitu¬ 
tional C appears to have been ejected out from the film at 
750 °C (based on the XRD measurements), which is unde¬ 
sirable. The reason for the ejection could be that this anneal¬ 
ing temperature is outside the window for substitutional C 
incorporation, 16 thus the metastable thin film tends to relax 
by precipitating substitutional C. 

V. SUMMARY 

We have studied the epitaxial growth of two groups of 
Ge-rich GeSiC films by low-temperature MBE. The films 
were extensively characterized using RHEED, TEM, Raman 
spectroscopy, x-ray diffraction, and spectroscopic ellipsom- 
etry. While some amount of C (up to 2-3 at %) is incorpo¬ 
rated in the GeSi alloys, excess C appears to cause 3D island 
formation and planar defects within the films. Increasing the 
amount of Si in the films enhances the tendency for planar 
defect formation, presumably due to lower mobility of Si 
atoms relative to Ge. 

Raman spectra of alloy films show that the effects of C 
on the strong Ge-Ge local modes are far less than the effects 
due to Si. We are unable to observe any systematic change in 
the Ge-Ge mode, whereas the Ge-Si mode appears to shift 
to lower frequency with the addition of smaller C. Spectro¬ 
scopic ellipsometry determinations of the film optical con¬ 
stants show that the primary effect of C is to reduce the 
strength of the E x critical point feature. 
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